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Hexagonal warping on optical conductivity of surface states in Topological Insulator 
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ARPES studies of the protected surface states in the Topological Insulator 512 Tes have revealed 
the existence of an important hexagonal warping term in its electronic band structure. This term 
distorts the shape of the Dirac cone from a circle at low energies to a snowflake shape at higher 
energies. We show that this implies important modifications of the interband optical transitions 
which no longer provide a constant universal background as seen in graphene. Rather the conduc- 
tivity shows a quasilinear increase with a slightly concave upward bending as energy is increased. 
Its slope increases with increasing magnitude of the hexagonal distortion as does the magnitude of 
the jump at the interband onset. The energy dependence of the density of states is also modified 
and deviates downward from linear with increasing energy. 

PACS numbers: 72.20.-i, 75.70.Tj, 78.67.-n 
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I. INTRODUCTION 

Topological Insulators are insulating in the bulk and 
have symmetry protected helical Dirac fermions on their 
surface with an odd number of Dirac points in the sur- 
face state Brillouin zoncii"^ The surface charge carriers 
are massless and relativistic with linear in momentum 
energy dispersion curves. Spin sensitive angular resolved 
photo emission spectroscopy (ARPES) also shows that 
their spins are locked to their momentum. The position 
of the chemical potential relative to the Dirac point and 
the gapped bulk bands is not as easily tuned as it is in 
graphenei^ but can be controlled by doping with Sn in 
(Bii-sSns)2Te^ or with Ca in Bi2-sCasTe^ with fur- 
ther dosing with NO2 molecules. The constant energy 
contours in Bi2Te^ as measured by ARPES are not cir- 
cles as they are in graphene but have an hexagonal dis- 
tortion which gives them snowflake shape. This geometry 
was modeled by Fuf^ with an unconventional hexagonal 
warping term in the bare band Hamiltonian of Bi2Te3 
with parameters fit to the measured Fermi surface. 

Optical spectroscopy has been a very powerful method 
to obtain valuable information on the charge dynamics 
of the Dirac fermions in graphene J^^— An experimental 
review was given by Orlita and PotemskiJ^ Usually it is 
the zero momentum q limit of the optical conductivity 
as a function of photon energy which is measured but 
very recently finite q's have also been measured using 
near field techniques li^""— Optics has also been used to 
study topological insulators i^"— In this paper we study 
how hexagonal warping in Bi2Te3 manifest in the optical 
conductivity, which we find is profoundly modified for the 
parameters determined in the work of Fu»ii 



II. FORMALISM 

The Hamiltonian used by Fi*ii to describe the surface 
states band structure near the F point in the surface Bril- 



louin zone is 

Ho ^ Vk{k^(Jy - kya^) + -(fc+ + kl)a^ + E^ik) (1) 

where Eoijz) = fi^k^ /{2m*) is a quadratic term which 
gives the Dirac fermionic dispersion curves an hour glass 
shape and provides particle-hole asymmetry. The Dirac 
fermion velocity to second order is Vk = vpi^ + ctk'^) with 
vf the usual Fermi velocity measured to be 2.55eV ■ A 
and a is a constant which is fit along with m* to the mea- 
sured band structure in the reference [ll|. The hexago- 
nal warping parameter A — 250eT^ ■ A^. The ax, cry, Gz 
are the Pauli matrices here referring to spin, while in 
graphene these would relate to pseudospin instead. Fi- 
nally fc± = fca; ± iky with kx , ky momentum along x and 
y axis respectively. The energy spectrum associated with 
the Hamiltonian [Eq. ([1])] is 



E±{k)=Eo{k)± 



v^k^ + X^k^cos^iW) (2) 



where 9 is the polar angle defining the direction of k in 
the two dimensional surface state Brillouin zone. The 
energy dispersion curves in Eq. ^ reduce to the well 
known linear law rkv^k of graphene when Eq is set zero 
along with A = 0, i.e. ignoring hexagonal warping. Since 
our primary interest here is getting a first understanding 
of how the warping term in Eq. ^ manifests itself in 
the dynamical conductivity of the surface helical Dirac 
fermions we will for simplicity from here on drop the 
Eoik) term which as we have said provides particle-hole 
asymmetry. 

In FIG. [T]we show a color plot for the constant energy 
contour associated with the dispersion curves [Eq. ([2])] 
as the energy is increased above that of the Dirac point 
the contour changes shape and shows greater hexagonal 
distortion displaying a snowfiake shape. The largest flake 
shown corresponds to a chemical potential of 250meV 
which is achieved when 5 = 0.67% in {Bii^sSns)2Te3. 
In graphene we would have circles for all energies and 



i 






1.5 



- 1 



).5 

) 

3.5 



0.5 

0.45 

0.4 

0.35 

0.3 

I ^ 0.25 
0.2 
0.15 
0.1 
0.05 




--1 



-2 -1.5 -1 -0.5 0.5 1 1.5 2 

K< r }K 



■1.5 

■2 




FIG. 1. (Color online) Constant energy contours for the dis- 
persion curves used to describe the bare bands in Bi2Tez 
with chemical potential ^ changed by doping with Sn in 
{Bii-.5Sns)2Te3 where 5 = .67% corresponds to a chemi- 
cal potential fj, — 250meV . The kx and ky axes are in the 
units of 0.1A~^ . Also shown is the surface state Brillouin zone 
identifying F, K and M points. 



in strained graphene we would have elliptical contours 
instead of circle. 

The Kubo formula for the xx component of the dy- 
namic conductivity a^xi^jj) as a function of photon energy 
w is given in terms of the matrix Matsubara Green's func- 
tion G{k,LUn) with w„ the Fermionic Matsubara imagi- 
nary frequency as 



e^ 1 
(Txxi^) = - — -r^ 



kdkde 



T ^ Tr{vxG{'k,uji)vxG(k,ujn + uji))i^^^^uj+iS (3) 



with e the charge on the electron, k the absolute value of 
the momentum with direction and kcut a cut off. Here T 
is the temperature with ujn = (2n + 1)ttT and w; = 21ttT 
the Fermion and Boson Matsubara frequencies ^21 n and 



I are integers and Tr is a trace. To get the conductivity 
which is a real frequency quantity, we needed to make an 
analytic continuation from imaginary iw„ to real to and 
S is infinitesimal. As written we have neglected vertex 
correction and so the factors Vx are simply the velocity 
components given by 

Vx — VkO'y + 3Afc^ cos(26')ct2 
Vy — —VkCTx — 3Afc sin(20)cr^ 



(4) 
(5) 



obtained directly from the Hamiltonian [Ec^. ([T])]. We 
have set all h factors equal to one. 

The matrix Green's function for the non-interacting 
bare band is given by 

G'o(k,iw„) 

1 

iuJn+ 11- Vk{kx(Ty - kyffx) - |(fc+ + kt)(Tz 



3, 



iun + A* + VkikxO'y — kyax) + \k'^ cos{39)a 
(iLUn + m)' - vlk^ - X^k^ cos2(30) 



(6) 



with 



(fc| + kty = {2ki - Qkxk'Y = 4fc'^ cos''(36') (7) 



It is convenient to rewrite the 6*0(14, ia;„) in terms of 
6*0(14, s, iuJn) defined as 



Go(k, s,itj„) 



1 



iwn + /i - s^Jvlk"^ + A2fc6 cos2(36') 

where s = ±1 and F^ defined as 

(— Wfefcsin6',Wfcfccos6', Afc'^cos(30)) 
''" V";^fc2 + A2fc6cos2(36l) 

This gives 

Go(k,ia;„) = :j ^(1 + sF^ ■ CT)Go(k, s,zw„) 

'^ s=± 



(8) 



(9) 



(10) 



III. SIMPLIFICATION OF EXPRESSION FOR 

We will be interested here only with the interband 
terms in which case the required trace gives 

Tr{vxG{\<i,uJi)vxG{\<i,i^i +cj„)} 

H{e) 



^fc2+A2fc6cOs2(36')] 



1 



1 



1 1 



(11) 



where 

H{e) == A2w2/c6[cos2(36I) -3cos6lcos(26')cos(36i) 

+9cos2(26l)]-|-w^fc2sin2 6' (12) 



But we know that 



^Ei 



1 



1 



ioji + fi - E- iuji + iuJn + fJ, - E+ 



1 



ME^)-f{E+) ^ f{E+)-f{E.) 



(13) 



and hence we obtain 



kdkdO- 



H{0) 



ivlk'^ + X'^k^cos'^iSe)] 

f{E_)-f{E+) , f{E+)-f{E^)^ 



JiU/V— >cj+'id 



ZWn — E+ + i/_ ZWn — £/_ + £/+ 

Here f{x) is the Fermi-Dirac distribution function given 
by f{x) = l/[exp(x/T — /i/T) + l], where we have ignored 
the Boltzman constant ks but will include it in the calcu- 
lation. We have verified that ayy{uj) = axx{'^)^ after an 
analytic continuation from imaginary to real Matsubara 
frequencies we obtain the final expression 

axx{^) = ^^-r^ I kdkdOTTTT-. — — x 



(15) 



47r2yo '"""^""Wik.e) 
f{E_)-f{E+) ^ f{E+)-f{E^) 
Lo-2^W{k,e)+i5 uj + 2^W{k,e)+i5 
where W{k, 9) = uf fc^ + A^fc^ cos2(36'). 



IV. ANALYTIC FORM FOR THE REAL PART 
OF CONDUCTIVITY 

The real part of the dynamic conductivity which is the 
absorptive part can be simplified further using the usual 
rule -^jTis ^ TJ ^ inS{uj). From Eq. (|15l) we get 

2 1 pkcut 

Reaxxiuj) = - — -r^ kdkd9H(e) x 

a; 47r^ Jq 



[f{E^)~f{E+)] 



[W{k,9)] 

which can be rewritten as 

e^ 1 /■'=■ 
Reaxxiuj) ^TTir d{k')de 

Zlo 47r Jo 



[-TT)5{uj-2y'W{k,e)) (16) 



w2/4 



[f{-^/2)-f{u:/2)]5{u:^2.JW{kJ)) (17) 

and we get 

e2[/(-^/2)-/(^/2)] 



Reaxxi^) 



247ra;2 
diP)H{ek^^) 



/o |A2fc6cos(30fe,„)sin0fc,^| '^^^^ 

where the thermal factor /{to) have been pulled out of 
the integral over momentum. The Dirac delta function 
has been used in the process. We wrote 



S{uj - 2y/W{k,e)) 



LuS{0 - Ok, 



\12X2k^cos{3ek,u>)sm0k,, 



(19) 



where 



1 , uj^/A-vlk\ 

±-arccos[±^ A^fc^' ] 



a;2/4-w?fc2 



±-{arccos[± 



±-{arccos[±^^^5^]+27r}. 



(20) 



This is our final expression for the absorptive part of the 
conductivity. Numerical results based on this expression 
are given in the next section. Before doing so however 
we present a similar expression for the density of states 
D{(jj) as a function of energy which could be measured in 
scanning tunneling microscopy (STM). By its definition 

D{uj)^ V/mrrG(k,iw„- > w + Z(5) (21) 

77- r J 



which can be reduced to 



D(u:) 
Im\ 



1 1 

IT 47r2 



kaut 



kdkde 



Lo- ^W{k,9)+i6 L0+ ^W{k,e)+i6' 
Taking the imaginary part gives 



(22) 



D{u) 



1 



kcut 



47r' JO 
+ (a; -^ -Lo) 

which works out to 

D{u) 



kdkdOSiuj - yJW{k,e))Q(uj) 



(23) 



kdk 



5^0M 



12^2^2;^ |fc6cos(30^^Jsin0^J 
+{uj -^ -w) (24) 

where Q{uj) is the Heaviside step function and we have 
used 



5{lj - ^w{k,e)) = 



^5{e~eij 



\■iX^k'^cos{W'^^J&lne'^ 



(25) 



where 6'J. ^ can be obtained from Eq. ([20|) by replacing lo 



with 2uj so we have 9',, 



7fc,2u 



NUMERICAL RESULTS 



In FIG. [21 we show our results for the real part of the 
optical conductivity Reaxx{^), in units of 2Tre'^/h, as a 
function of photon energy uj for a Bi2Te3 doped with 
Sn at level S = .67% (see references (ja) and (jllI )) which 
corresponds to a chemical potential fj, = 250meV and 
all other parameters determined in the fit by Fu.— We 
show four values of /x. In all cases the threshold for the 
start of the interband transitions is sharp and occurs at 
Lo = 2/x, as it would in graphene. The missing optical 
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FIG. 2. (Color online) The real part of the optical conduc- 
tivity axx{io) as a function of photon energy lj in meV (in 
units of 2-Ke^ /h) for a case which corresponds approximately 
to S — 0.67% Sn doping with chemical potential fi = 0.25ey. 
We show 4 values of fi. In all cases finite /i transfers optical 
spectral weight from the interband transitions to the intra- 
band. These, not shown here, provide a Drude like contribu- 
tion at oj near zero. 



spectral weight in the interband transition is accompa- 
nied with an increase in the intraband (Drude) optical 
spectral weight. This is not shown in our picture. Be- 
cause we have not included any scattering processes in 
our work, the Drude manifests as a Dirac delta function 
at w = and does not overlap with the interband contri- 
bution which we emphasize here. At small values of uj, 
the value of Rea(u)) is rather flat and takes on precisely 
the value expected for graphene without the degeneracy 
factor 5 = 4, which counts spin and valley degrees of free- 
dom. For a topological insulator there is only one Dirac 
cone and spin is no longer degenerate. We also note that 
the background value is independent of material param- 
eters such as the Fermi velocity. But this is no longer 
the case for a topological insulator. As to is increased 
whatever the value of ^ the conductivity increases rather 
rapidly above its universal background value and shows 
concave upward behavior. This is traced to the changes 
in fermi velocity of Eq. (|4|) and Eq. ([5|) due to the warping 
term proportional to A and to the change in quasiparticle 
band structure. In FIG. [3] we show how Rea{uj) v.s. uj 
is changed as A is changed. Here and also in FIG. |4]the 
A has been multiplied by the cube of the typical Fermi 
momentum, which is 0.1A~^. So A = 0.2 here would cor- 
respond to A = 200eF • A^. For reference the case A = 
given as the solid black curve, corresponds precisely to 
graphene except for value of the degeneracy factor g. In 
this case the universal background, well known in the 
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FIG. 3. (Color online) The real part of the optical conductiv- 
ity CTa:i,(cj) in units of 2-Ke^ /h as a function of photon energy 
for 4 values of the magnitude of the hexagonal warping term 
(A) in the Hamiltonian (1). The solid black curve is for com- 
parison. In this case A = and the model reduces to the 
contribution to the optics of a single spin and single valley 
Dirac cone of graphene. In the inset we show the increase in 
the value of the jump at the interband absorption edge with 
increasing A. 



graphene literature, is recovered. As /i is increased the 
optical spectral weight lost in the background is trans- 
ferred to the intraband contribution not shown above. 
Increasing A changes the band structure and no optical 
sum rule applies. As A is increased at fixed value of 
chemical potential fj, the magnitude of the interband on- 
set {cT onset) (which remains at 2^) increases as shown in 
the inset. In addition its increase at oj > 2/x becomes 
even more rapid. While in the black dashed curve it is 
reasonably linear, the dotted curve for A = 0.2 has ac- 
quired a significant upward curvature. These deviations 
from the universal background of graphene is the sig- 
nature in optics of the hexagonal warping term in our 
Hamiltonian [Eq. ([T])]. We have checked and found that 
these curves do not scale onto each other. The increase 
in Rea[(jj) with w above the universal background value 
is to be contrasted with the behavior of the quasipar- 
ticle density of states D{(jj) given by Eq. (24). As in 
graphene, the case A = gives D{lu) proportional to |w|. 
As A is increased however D{uj) starts to deviate from lin- 
earity and, as we see in FIG. |31 is progressively reduced 



VI. SUMMARY AND CONCLUSION 
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FIG. 4. (Color online) The density of state D{u)) as a func- 
tion of u). We show these for values of the hexagonal distor- 
tion term A = (solid black line), A = 0.1 (dashed red line) 
and A = 0.2 (dotted blue line). The solid black straight line 
corresponds to graphene without spin and valley degeneracy. 
The Fermi velocity vf was taken to be that appropriate to 
Bi2Te3. The vertical dashed black line indicates a chemical 
potential fi = 0.25ey. The left inset gives the value of /i as 
a function of A for fixed number of charge carriers which is 
around l.l*10^®/m^. The right inset shows the Fermi surface 
at fixed value of /i for A = (black), A = 0.1 (red), A = 0.2 
(blue). With no hexagonal distortion it is circular and distorts 
to a snowfiake as A increases. 



below the solid black curve. This is easily understood 
with the help of the right inset where we plot the con- 
stant energy contours for w = /i = 0.25eF for the three 
value of A considered. For A = we get the black circle 
of graphene theory. As A increases this contour distorts 
into a snowfiake pattern (blue curve) which is however 
completely contained inside the black circle. Of course, 
to keep the number of charge carriers the same we need 
to increase the chemical potential with increasing value 
of the warping parameter A as we show in the left inset. 
What is plotted is the value of A at fixed value of the 
number of charge carriers, which is around 1.1 * 10^^ /m^. 



Helical Dirac fermions exist at the surface of a topo- 
logical insulator (TI). These charge carriers have some 
similarity and also differences with the well known chiral 
Dirac fermions in graphene. An important difference is a 
degeneracy factor g of four which comes from the valley 
and spin degrees of freedom of graphene not applicable 
in TI. Another important difference, well investigated in 
the case of Bi2Te^ doped with Sn, is the hexagonal dis- 
tortion seen in ARPES. Here we have studied how such 
a term changes optical properties. For realistic values 
of the warping parameter we found large changes in the 
interband transitions. 

A third difference is that, graphene involves pseudospin 
related to the sublattice degeneracy in its two atoms per 
unit cell crystal structure rather than real spin. Further- 
more in graphene the bands are spin degenerate, while 
in a topological insulator momentum and spin are locked 
with x-y component of real spin oriented perpendicular 
to its 2-D momentum k, with clockwise and anticlockwise 
orientation in conduction and valence band respectively. 

The universal flat background observed in graphenei^ 
remains at small photon energies although modified by a 
factor of 4 because the valley spin degeneracy no longer 
applies. As lo increases large modifications in the ef- 
fect of the interband transitions on the conductivity are 
noted, and these encode the information on the hexago- 
nal warping of the Dirac cone cross-section leading to a 
snowfiake pattern. Instead of being flat Rea{ui) increases 
in a quasilinear fashion with a concave upward bent. The 
magnitude of this linear increase becomes larger with the 
magnitude of the hexagonal warping term as does the 
value of the jump in the conductivity at the threshold 
of twice the chemical potential (w = 2/x). At the same 
time, we find that the density of state remains linear only 
at small w and starts to fall below this linear behavior 
at the energy where the conductivity also starts to show 
its deviation from a constant background value. While 
the conductivity curves are bent upward due to fermi ve- 
locity features, the density of state bends downward a 
prediction that could be verified in combined optics and 
scanning tunneling spectroscopy(STS) experiment. 
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